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Abstract 


Recent technological advances in airplane construction techniques and 
materials employing bonded and milled aluminum skins and composite materials 
allow for the production of aerodynamic surfaces without significant wav i ness 
and roughness, permitting long runs of natural laminar flow (NLF). These 
advances lead to excellent opportunities for reducing the drag of aircraft by 
increasing the extent of NLF. The present research effort seeks to refine and 
validate computational design tools for use in the design of axisymmetric and 
nonaxi symmetric natural -laminar-flow bodies. The principal tasks of the 
investigation involve fuselage body shaping using a computational design proce- 
dure. 

Under Phase I SBIR funding for this research, analytical methods were refined 
and exploratory calculations were conducted to predict laminar boundary-layer 
behavior on selected body shapes. Using a low-order surface-singularity aero- 
dynamic analysis program pressure distribution, boundary-layer development, 
transition location and drag coefficient have been obtained for a number of 
body shapes including a representative business-aircraft fuselage. Extensive 
runs of laminar flow were predicted in regions of favorable pressure gradient 
on smooth body surfaces. A computational design procedure was developed to 
obtain a body shape with minimum drag coefficient having large extent of NLF. 

Some preliminary results from the design efforts have been obtained and further 
work is underway. 

The proposed study has widespread commercial applications. A significant 
reduction in the drag produced by any airplane can be obtained when extensive 
runs of natural laminar flow are achieved on its fuselage resulting in improved 
airplane performance and efficiency. 



Nomenclature 


Drag coefficient 
Pressure coefficient 
Maximum diameter 

Fineness ratio, length/maximum diameter 

Boundary-layer shape factor, 6*/0 

Nondimensional curvature at x m 

Body length 

Mach number 

Natural Laminar Flow 

Logarithmic exponent of T-S wave growth ratio 

Profile radius at x n - 

Radius of curvature at the nose 

Reynolds number based on free-stream conditions and body 
length 

Reynolds number based on local conditions and surface 
length 

Reynolds number based on local conditions and axial length 

Body volume Reynolds number based on free-stream conditions 
and V 
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Reynolds number based on local conditions and boundary-layer 
momentum thickness 

Nondimensional profile radius of x^, 2f r R^/L 
Nondimensional radius of curvature at the nose, 4x |T| f r R r| /L 
Nondimensional profile radius at x 
Profile slope at x^ 

Surface length, starting at the nose 
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Free-stream velocity 
Local velocity 
Body volume 

Axial coordinate, starting at the nose 
Axial location of inflection point 
Axial location of maximum diameter 

Nondimensional axial location of maximum diameter, XVL 

m 

Nondimensional axial location of inflection point, X^/L 

Nondimensional axial length coordinate of boundary-layer 
separation 

Nondimensional axial location, X/L 
Nondimensional thickness distribution 
Angle of attack (degrees) 

Boundary-layer displacement thickness 
Boundary- layer momentum thickness 

Half trail ing-edge angle at the tail (see fig. 15) in degrees 
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1. Introduction 


In recent years, airplane construction material and fabrication methods 
have improved greatly, resulting in the production of airframe surfaces which 
accurately match the design shape. Recent flight tests (refs. 1 and 2) have 
demonstrated that extensive runs of laminar boundary layer flow can be obtained 
over regions of favorable pressure gradient on smooth airplane surfaces and 
provide a significant reduction in profile drag. A major portion of the past 
research effort for achieving and maintaining natural laminar flow (NLF) has 
been focused on aircraft wings. Fuselage shaping to increase the extent of 
NLF has received much less attention in the literature except for sailplane 
bodies and hydrodynamic bodies. 

Althaus (ref. 3) conducted experimental investigations in order to show 
the possible reduction in drag of sailplane fuselages and study the flow inter- 
action of the fuselage with a wing and the influence of various wing positions. 
In reference 4 an analytical and experimental study on NLF nacelles demonstrates 
the presence of significant regions of laminar flow on the surface of a turbo- 
fan engine nacelle. Sub-scale wind-tunnel testing of the NLF nacelle indicated 
a potential of 1. 5-2.0 percent reduction in total airplane drag by maintaining 
laminar boundary- layer flow over extensive regions of the external nacelle 
surface. The pay-off in terms of airplane profile drag reduction can be much 
larger when considering the application of low-drag NLF design to airplane 
fuselages. The importance of fuselage skin-friction drag is clearly indicated 
in figure 1, in which a profile drag breakdown is shown for a typical transport 
jet (ref. 5). The fuselage generates less than 50 percent of the profile drag 
for the all-turbulent airplane. However, it is estimated that the contribution 
of the fuselage to airplane profile drag increases to more than 70 percent of 
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the total profile drag if extensive regions of natural laminar flow are 
achieved on the wing and tail surfaces. The present study will investigate 
the design of fuselage shapes that result in extensive regions of NLF at 
conditions corresponding to typical cruise Reynolds numbers encountered by 
touring, business and transport airplanes. 

Carmichael (ref. 6) did an experimental study on a body of fineness 
ratio of 3.33 developed by revolving the coordinates of a NACA-66 laminar flow 
airfoil about the longitudinal axis. This tailboomed body (called the Dolphin) 
was tested over a Reynolds number range of 20 million to 30 million (based on 
body length) and transition length Reynolds numbers of 14 million to 18 million 
were obtained. Boundary-layer transition from laminar to turbulent flow 
apparently occured beyond the point of maximum thickness. These results 
indicate that a low fineness ratio and a proper shape can produce a strong 
favorable pressure graident (and therefore a strong flow acceleration) on the 
forebody of the configuration. As a result the boundary layer stays laminar 
over an appreciable distance. 

Dalton and Zedan (ref. 7) presented an inverse method to design low-drag 
axisymmetric body shapes. A prescribed surface-velocity distribution is input 
and the corresponding body shape is computed. The method is based on represent- 
ing the body of revolution by a source distribution of variable intensity along 
its axis. Dalton and Zedan applied the method to carry out a design study for 
bodies at a volume Reynolds number R y = 50 million. This Reynolds number is 
representative for large torpedo and small airship applications and it 
corresponds to a length Reynolds number R L of approximately 200 million 
depending on body shape. 

Forebody shapes of missiles designed specifically for long runs of NLF at 
compressible free-stream velocities have been studied in reference 8. In 
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reference 8, results computed for three forebody shapes of different fineness 
ratios (1.0, 1.5 and 2.0) indicate that laminar flow can be obtained for high 
Reynolds numbers and high subsonic Mach numbers. The design condition for 
these relatively blunt forebody shapes was a unit Reynolds number of 40 million 
per foot and a Mach number of 0.75. However, the Reynolds number range which 
is encountered by most airplanes is in the range of 1-3 million per foot. 
Therefore, the results of reference 8 are not very useful for the design of 
fuselage forebodies for subsonic airplanes. 

Parsons and his coworkers (refs. 9 and 10) have given a procedure to 
design axisymmetric bodies for minimum drag for hydrodynamic applications. 

Drag reduction is accomplished solely through manipulation of the vehicle shape. 
The optimization problem is formulated as a nongradient search in a finite 
constrained parameter space. Two classes of bodies, described by five and 
eight parameters, are considered. The requirement for nonseparating flow 
represents an additional constraint on the optimization problem. The axisym- 
metric bodies are represented by axial singularity distribution and the drag 
coefficient is evaluated using Young's formula. The results show that signifi- 
cant drag reduction is possible through shape manipulation. Pinebrook (refs. 

11 and 12) devised a technique, based on an evolution strategy, to minimize 
the drag of an axisymmetric body with a given maximum body diameter and fineness 
ratio. The body profile is described by continuous first order axial 
singularity distribution defined at 21 points. The gradual body profile 
changes are effected through a process derived from the evolution strategy. The 
drag is calculated from the momentum deficit in the boundary layer at the 
end of the body using Young's formula (ref. 13). The bodies designed by 
Pinebrook had laminar flow onlyupto 3 percent of the body length from the nose 
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at which location the boundary layer was tripped. 

According to von Karman (ref. 14) not every axi symmetric body can be 
represented accurately by an axial source distribution. Also, recently Hess 
(ref. 15) pointed out that the line sources cannot extend to the ends of the 
body at finite strengths otherwise the velocities there will be infinite. These 
problems with the axial singularities are eliminated if one uses panel methods 
with which the body is represented by aerodynamic surface panels. Here, a 
computational procedure was developed to obtain fuselage geometries for consider- 
able extent of laminar flow and hence for low skin-friction drag. 

Computational Design Procedure 

The computational design procedure used to obtain "natural laminar flow 
fuselage" geometries is described in the flow chart (see fig. 2). Initial 
values of the design variables describing the body shape are input along with 
the Reynolds number (based on length), the nondimensional length and the 
fineness ratio. The axi symmetric body is described by seven parameters. The 
expressions to obtain the body shape are given in the Appendix. 

A number of constraints are imposed on the design parameters in order to 
generate designs which are realistic and practical. The geometric constraints 
are given below. 

°^n 

2 . 0 < k x 

3. 0 < x < x- < 1 

m l 

4. 0 < r- £ 1 

5. 0 £ s i 

6. 5° < $ < 80° 

7. No inflection on forebody, midbody and afterbody except at x^. 
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These conditions are taken care of by choosing proper upper and lower bounds 
for the design variables. 

The other constraint is that the separation takes place very near the 
trailing edge. The following constraint is imposed in the optimization proce- 
dure: 

0.95 < x 
= sep 

The objective function is taken to be body profile drag coefficient which is 
obtained by Young's formula (ref. 13). This objective function is to be 
minimized subjected to the constraints which are given above. The optimizer 
computes the gradients of the objective function and then, using either a 
conjugate direction method or a method of feasible direction, determines a 
linear search direction, along which a new constrained variable is constructed. 

An improved or minimum feasible objective functional value is evaluated 
and a series of proposed updated design variables are calculated. The objective 
function and the constrained function are evaluated using the updated design 
variables, interpolating over the range of feasible proposed design variables 
resulting in a minimum value of the objective function. The results are tested 
against a convergence criteria. The procedure will stop if the convergence 
criteria is satisfied giving a body shape with minimum drag or maximum transi- 
tion length satisfying the separation constraint. If the convergence criteria 
is not satisfied the design parameters go through the analyzer again resulting 
in a new set of design variables and the procedure is repeated until a final 
shape is obtained. The method involves a constrained minimization procedure 
(ref. 25) coupled with an aerodynamic analysis program based on a low-order 
surface singularity method named "VSAERO" (ref. 16). Pressure distribution and 
velocity distribution can be computed by the aerodynamic analysis program which 
uses surface singularity panels to represent the body. The program also calcu- 
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lates the effect of the viscous boundary layer adjacent to the body. Integral 
methods are used to predict the boundary-layer development. The laminar part 
of the boundary layer is calculated by Thwaite's method with Curie's modifica- 
tion. Boundary-layer transition is predicted by Granville's procedure. Nash 
and Hick's method is used for the turbulent boundary-layer calculations. 
Laminar separation/turbulent reattachment calculations are done empirically 
using Gaster's measurements. 


Aerodynamic Analysis 

To validate the surface singularity method, inviscid pressure distributions 
have been obtained for the following configurations: 

(a) a sailplane body, 

(b) an axi symmetri c body with an ogival nose, a cylindrical 
center body and a flared afterbody, 

(c) an ellipsoid of revolution, 

(d) a body of revolution with a long favorable pressure 
gradient, 

(e) a low-drag body of revolution considered by Parsons 
et al . (ref. 10) 

(f) a representative business aircraft fuselage, and 

(g) a body of revolution whose maximum diameter and length 
correspond to those of the configuration (f). 

In reference 3, Althaus presents wind-tunnel measurements on bodies and 
wing-body combinations obtained at the University of Stuttgart. The measurements 
include lift and drag characteristics, surface velocity distributions, and 
boundary- layer transition locations. Each of the bodies of revolution had a fineness 
ratio of 10 and a length of 6.56 ft (2.0 m) and the measurements were obtained at 
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a Reynolds number (based on body length) = 7.1 million. Althaus' body 
shape 2 was modeled by surface singularity panels to obtain velocity distribu- 
tions. In figure 3, the theoretical velocity distributions are compared with 
surface velocity measurements. In the figure, the body shape is also shown 
for this fuselage model. In the pressure recovery region theoretical results 
do not agree very well with the wind-tunnel data. Also, at X/L = 0.12 a kink 
appears in the calculated velocity distribution. This disturbance is produced 
when the body coordinates as given by Althaus in reference 3 are used to model 
the body shape. Consequently, the prediction of the boundary- layer transition 
location does not agree very well with the experimental result. In the tunnel, 
transition was measured at X/L = 0.36 at R^ = 7.1 million. The Granville 
criterion predicted transition location at X/L = 0.29 which is 
slightly forward compared to the experimental result. This can be expected 
because of the velocity discrepancies in the pressure-recovery region. These 
discrepancies appear to be caused by wind-tunnel blockage effects, which 
produce higher surface velocities near the location of maximum body diameter. 

Very good agreement between measured and calculated surface pressure 
coefficients was obtained for a body of revolution (configuration b) tested by 
Fox in the NASA Langley high-speed 7- by 10-foot wind tunnel (ref. 17). In 
figure 4, the results are shown for Fox's configuration 5, an axisymmetric body 
composed of an ogival nose, a cylindrical centerbody and a flared afterbody. 
Surface pressure data are plotted as function of the orifice location for the 
orifices on the longitudinal meridian. By using the surface singularity 
method, the pressure coefficients are obtained and they are also presented in 
figure 4. The method allows modeling of the wake which originates from the 
edge of the blunt base. However, for simplicity when the blunt base and the 
sting were replaced with a 10-degree cone (the flared afterbody has an identical 
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slope of 10 degrees) the comparison of theoretical results with experimental 
data showed very good agreement. 

Measurements of pressure and boundary-layer transition were made on the 
ellipsoid of revolution (configuration c) in wind tunnels and boundary-layer 
transition was also measured in flight (ref. 18). Figure 5 presents the surface 
pressure distribution and transition location for the ellipsoid of revolution of 
fineness ratio 9. The viscous calculations were done at R^ = 13.98 million. 

As shown in the figure, although the pressure distributions compare well, the 
computed transition location does not agree very well with the experimental data. 
The transition length Reynolds number was 5.08 million. Granville's boundary- 
layer transition criterion, which is used in the present computation, predicts 
transition further downstream than was measured in the experiment (ref. 

18). The same behavior occured when the Reynolds number was increased to 22.03 
million. Kaups (ref. 19) also computed transition location for this ellipsoid 
by using Granville's criterion and also arrived at the fact that at these two 
Reynolds numbers, transition is predicted at points further downstream than at 
the locations observed in the experiments (in flight and wind tunnel). Kaups 

Q 

concluded that none of the transition prediction methods (including the e -method 
of Smith and Gamberoni) gave consistent satisfactory answers for the ellipsoid. 

It appears that for bodies with flat pressure distributions considerable amount 
of uncertainty exists in the computation of transition location. t 

Next, a low-drag body of revolution with fineness ratio 4.5, with a long 
favorable pressure gradient forebody (configuration d) was analyzed. Hansen 
and Hoyt (ref. 20) studied this body experimentally and measured drag and 
intermittency and calculated the surface-pressure distribution. In figure 

tin private communication. Dr. Pfenninger has indicated that the F94 boundary- 
layer transition data appears to be adversely affected by excessive engine sound 
levels. 
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6, comparison is made between Hansen and Hoyt's results and the predictions 
obtained by the present method. In the present computations transition is 
predicted at X/L = 0.67 due to laminar separation for = 10.86 million. 

Hansen and Hoyt predicted transition at X/L = 0.68. 

Another low-drag body called X-35 (configuration e) studied by Parsons 
et al . (ref. 10) was also analyzed. The predicted surface-velocity distribution 
agrees very well with that presented in reference 10 (see fig. 7). Parsons 
uses the Michel-e correlation to predict boundary-layer transition. The 

9 

Michel-e correlation provides a relation between R Q and R at transition 

U S 

(ref. 21): 

R 0 = 1.174 (1 + 22400/R s )R s 0 * 46 (1) 

The results in figure 8 indicate that no catastrophic Tollmien-Schl ichting 
growth is predicted and transition due to laminar separation occurs at X/L = 

0.68. Granville's transition criterion is more conservative for the X-35 body 
with transition predicted at X/L = 0.35 for R L = 37.14 million. This discrepancy 
in the prediction of transition location has been further analyzed by applying 
the H-R x boundary-layer transition criterion by Wazzan, Gazley and Smith (ref. 

22) and also by performing a linear stability analysis of the laminar boundary 
layer. The H-R x method correlates the boundary- layer shape factor H = 6*/6 and 
R s at transition as follows: 

log[R $ (e 9 )] = -40.4557 + 64.8066 H - 26.7538 H 2 

3 (2) 

+ 3.3819 H J for 2.1 < H < 2.8 

The results in figure 9 show that according to this method boundary-layer 
transition occurs at X/L = 0.25. The linear stability analysis of the laminar 
boundary layer has been performed using the SALLY code (ref. 23). Input to this 
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program is provided by a modified version of the Harris finite-difference 
boundary- layer code (ref. 24). In figure 10, the logarithmic disturbance 
amplitude ratio or "n-factor" is plotted as a function of the nondimensional 
axial distance X/L for a range of Tollmien-Schlichting disturbance frequencies. 

The envelope of these curves shows that an n-factor of 9 is reached at X/L = 0.185, 
while at X/L = 0.25 the n-factor is about 12.5. In summary, it appears that the 
Michel-e 9 criterion as used by Parsons et al . (ref. 10) provides transition 
results which are too optimistic. 

Next, a representative business aircraft fuselage of fineness ratio about 
6 (configuration f)was analyzed. Inviscid pressure distributions on the upper 
and the lower surface were calculated and are presented in figure 11. Laminar- 
turbulent transition prediction was made using Granville's criterion at a unit 
Reynolds number of 1 million per foot (R L = 40.86 million) and it is indicated 
in the figure (transition Reynolds number range is 4.1 million - 16.0 million). 

It is to be noted that it is not possible to analyze these kinds of practical 
nonaxi symmetri c bodies by using axial singularity distributions. Furthermore, 
the transition prediction assumed no three-dimensional boundary- layer stability 
effects. While this assumption might limit the validity of the present transition 
prediction, the value of that limit (say in terms of Reynolds number) is not 
known. The validity of this prediction remains to be checked with experimental data. 

Finally a body of revolution whose maximum diameter and length (configuration 
g) correspond to those of the representative business aircraft fuselage config- 
uration (configuration f) is considered. Also this configuration can be 
thought of as a special case obtained by setting the longitudinal camber equal 
to zero in the more generalized configuration given by configuration f. Predicted 
inviscid pressure distributions and transition location are presented in figure 
12 for zero angle of attack. 
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Results from the Design Procedure 

Preliminary results obtained by the optimization procedure are discussed 
through an example and are given in Table 1. The analysis is for an angle of attack 
of zero and the flow is assumed to be incompressible. Initially the 
axisymmetric body is modeled by a set of initial values of the design para- 
meters. The fineness ratio is fixed at 6.14 and the design Reynolds number 
(based on length) R^ is 40.86 million. The upper and the lower bounds of the 
design variables are also input. In the present calculation profile drag 
coefficient is chosen as the objective function to be minimized. The design 
parameters at the end of each iteration are presented in Table 1 along with the 
objective function. For the above example the design procedure converged at 
the end of the 8th iteration. Since the initial values of the design variables 
were not too far from the optimized design variables the process converged in 
8 iterations. Judicious choice of the starting solution helps in achieving the 
final optimized solution in a lesser number of iterations. The initial body 
shape, given by the design parameters in the first row of Table 1, is sketched 
in figure 13. The pressure distribution for the body is also shown in the 
figure. For this initial body the transition location was computed at X/L = 

0.36. The improvements in the transition location through the next few 
iterations along with the final converged value are indicated in Table 1. We 
see from Table 1 that the transition location is pushed further aft for the 
final body shape and accordingly the drag is reduced. The final body shape 
along with the pressure distribution are sketched in figure 14. 

As explained earlier for computing the transition locations Granville's 
criterion is used in the analysis. Since to date there exists no satisfactory 
method to predict the transitional region on an axisymmetric body Granville's 
criterion has been used in the optimization cycle. Since one can predict 
growth of Tollmien-Schl ichting waves downstream of the neutral point in the 
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laminar boundary layer using linear stability analysis it is logical to include 
it in the design process itself. But in view of the computational costs and 
complexities, direct incorporation of the boundary-layer stability analysis in 
the optimization procedurewas avoided and Granville's criterion was used for 
predicting the transition. 


Summary 

Present airplane construction techniques result in the production of smooth 
and accurate aerodynamic surfaces over which long runs of natural laminar 
boundary- layer flow (NLF) can be obtained. A major portion of the past research 
effort for achieving NLF has been focused on airplane lifting surfaces. However, 
fuselage shaping to achieve considerable regions of NLF has received limited 
attention. Recent introduction of business and commuter airplanes with low- 
drag NLF lifting surfaces justifies study of the feasibility of NLF over fuselage 
surfaces and thereby providing possibilities to decrease airplane profile drag. 

In Phase I of the research effort, a low-order surface-singularity analysis 
was used to obtain surface-pressure distribution over a selected number of 
body shapes, including a representative nonaxi symmetric business-aircraft 
fuselage. Using an integral boundary- layer method, predictions of transition 
location and drag coefficient have been obtained for these body shapes. Exten- 
sive runs of NLF are predicted over regions of favorable pressure gradient on 
smooth surfaces of both axisymmetric and nonaxi symmetric bodies at low angles 
of attack. A computational design procedure coupled with an aerodynamic analysis 
program was developed to obtain a body shape with large extents of NLF and 
minimum drag. 

Comparisons of boundary- layer transition location predictions with limited 
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available experimental results for several bodies show poor correlation between 
calculated and measured transition locations. Therefore, it is proposed to use 
linear boundary-layer stability theory in the analysis and design of axi symmetric 
bodies with favorable pressure gradients for considerable runs of NLF. Presently, 
analytical tools are not available to calculate three-dimensional boundary- layer 
velocity profiles and transition location accurately under nonaxi symmetric flow 
conditions. Also, very limited high quality boundary-layer data under axisym- 
metric and nonaxi symmetric flow conditions are available at high Reynolds 
numbers to validate future developments in boundary- layer (stability) theory. 

Based upon the results of Phase I it is recommended that a large-scale 
high-Reynolds-number wind-tunnel experiment be conducted to obtain high-quality 
boundary- layer velocity and transition data on bodies of revolution at various 
angles of attack. The wind-tunnel experiment which would be designed using the 
optimization method developed in Phase I and linear boundary-layer stability 
theory, would provide needed experimental data for theory validation and also 
demonstrate the achievabil ity of extensive runs of NLF on commuter airplane and 
business- jet type of fuselages. 
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Appendix 


The seven design parameters describing the body shape are taken to be: 

x = nondimensional axial location of maximum diameter D 
m 

defined as X/L, 
m 

kj = nondimensional curvature at X m defined as (-2x m f r ) Kj L, 

r = nondimensional radius of curvature at nose defined as 
n 

^ 4x m W L * 

r. = nondimensional profile radius at X^ defined as 2f r R^/L, 

s^ = nondimensional profile slope at X^ defined as 

[- 2 V*i - 

x^ = nondimensional axial location of inflection point 

defined as X^/L, 

<f > = half trail ing-edge angle of the profile, 

and X^, Kj, R p , R^, and X^ are respectively the axial location of the 
maximum body diameter, curvature at X , radius of curvature at the nose, radius 
at X^, slope at X^ and axial location of the inflection point. 

The body is devided into 3 sections each described by a low degree 

polynomial. The complete body profile is continuous through the second 

derivative to avoid local regions of highly accelerated flow. The expressions 

representing the forebody and the midbody are taken from reference 10. The 

expressions representing the afterbody were derived according to the present 
needs. The body has a rounded nose and a pointed conical afterbody. 

The forebody (0 £ X £ X ) is described by a fourth-degree polynomial, the 
midbody (X < X _< X.) by a fifth degree polynomial and the afterbody by a fourth 
degree polynomial. The profile has a finite slope at X = L allowing a finite 
trail ing-edge angle. The dimensional parameters are shown in the figure 15. 

The analytical expressions for the nondimensional radius of the axisymmetric 
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body are given by 


r(x) = (^-)[r n F j(x) + kl F 2 (x) + G(x )] 1/2 

for 0 £ X < X m (forebody) 
where x = X/X 

m 

F^x) = - 2 x(x-l ) 3 
F 2 (x) = -x 2 (x - l ) 2 
G(x) = x 2 (3x 2 - 8 x + 6 ) 

For ( m i dbody ) 


r(x) = (^p) {r. + (1 - r.JCkj Fj(x) + Si F 2 (x) + G(x)]} 
r m ^ 

X, - X 

where x = ( x _ ) 

i m 

Fj (x) = - | x 3 (x - l) 2 
F 2 (x) = x - x 3 (3x 2 - 8x + 6) 


G(x) = x 3 ( 6 x 2 - 15 x + 10 ) 

ki 


k i = [$■> " 1] 

i m x m 
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For X.j jc X L (afterbody) 

r (x) = s. x ( 1 - x 3 ) - s. x 2 (2x - 3 ) (x - 1 ) 
L 'a 

+ x 2 ( 3 x 2 - 8 x + 6 ) 

where 

, . (1 - r j) - v s i 

f a < x i ' x m> r i 
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, 2f r (x i - x m> tan * , 
= ' [T^TJ 1 


A14 


f r = L/D 
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Convergence of the Design Variables (SPKB1 
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Reference area for C n calculation was frontal area. 
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Figure 1. Build-Up of Profile Drag for All-Turbulent Airbus A300 and A300 with Laminar Lifting 
Surfaces. 
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Figure 2. Flow Chart of Design Procedure for Minimum Drag Natural Laminar Flow Fuselage. 
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Figure 3. Velocity Distribution, Geometry and Transition Location for 

Fuselage Shape 2; a = 0°, R. = 7.1 million and Incompressible Flow. 
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Figure 4. Surface Pressure Distribution and Geometry for Model 
Configuration 5 at M = 0.404 and a = 0.06°. 




24 


Figure 5. Surface Pressure Distribution, Geometry and Transition Location 
for an Ellipsoid of Revolution of Fineness Ratio 9; a = 0°, 

R. = 13.98 million and Incompressible Flow. 
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Figure 6. Pressure Distribution, Geometry and Transition Location 

for a Body of Revolution with Fineness Ratio 4.5 ; a = 0°, R^ = 10.48 

million and Incompressible Flow. 
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Fiqure 7 Surface Velocity Distribution, Geometry and Transition Location for the X-35 Body 
of Revolution, a = 0°, R, = 37.14 million and Incompressible Flow. 
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Figure 9. Boundary-Layer Development and Predicted Boundary-Layer Transition for the 

X-35 Body at a = 0°, R, = 37.14 million using H-R x Method, Incompressible Flow. 
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Figure 10. Logarithmic Amplification of Several Tollmien Schlichting Disturbance 

Frequencies for the X-35 Body of Revolution at a = 0°, R. = 37.14 million 
and Incompressible Flow. L 
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Figure 11. Surface Pressure Distribution, Geometry and Transition Location 
J on a Representative Business Aircraft Fuselage; a = 0°, 

r. = 40.86 million and Incompressible Flow. 
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Figure 12. Surface Pressure Distribution, Geometry and Transition Location for a 

Body of Revolution in the Fineness Ratio 6.14; a = 0°, R, = 40.86 million 
and Incompressible Flow. L 




Figure 13. Surface Pressure Distribution, Geometry and Transition Location 
for Initial Body of Revolution, Fineness Ratio 6.14 and 
R = 40.86 million and Incompressible Flow. 
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Figure 14. Surface Pressure Distribution, Geometry and Transition Location for 
Optimum Body of Revolution, Fineness Ratio 6.14, R, = 40.86 million 
and Incompressible Flow. L 
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Figure 15 
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